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It is shown that the multiplicity distribution associated to self-shadowing events 
satisfies an universal equation in terms of the minimum bias distribution. The 
number of elementary collisions, basic in the structure of any multiple scattering, 
is reduced if a collective effect like percolation of strings takes place. The main 
consequences of the percolation are briefly discussed. 



1 Introduction 

Most of the exciting and forthcoming data on particle production on hadron- 
hadron and nucleus- nucleus collisions can be understood by conventional physics. 
This does not mean lack of excitement. Indeed, in this paper we show first as 
an example a surprising general law on multiplicity distributions associated to 
events which are shadowed by themselves, obtained using nothing but unitarity 

Also, some data of the forthcoming experiments in the Relativistic Heavy 
Ion Collider (RHIC) and in the Large Hadron Collider (LHC) can give us many 
surprises related to new physics. As an example of this, some effects of the 
percolation of stringscra, that can be detected already at RHIC, are discussed 
in this paper. 



2 Multiplicity Distributions Associated to Events Shadowed by 
Themselves. 

A multiple collision can be seen as a superposition of many elementary colli- 
sions, and it is worth to ask for the propagation of properties existing at the 
level of elementary collisions. In this way, for example, it is well known the 
A— behavior of the hadron-nucleus cross sections for events which are shadowed 
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by themselveslffl. Indeed, for these events, denoted by C-events, the resulting 
h — A event is of type C only and only if at least one of the elementary collisions 
if of type C. It is said, the algebra 



C x no C = C 
no C x no C = no C 
CxC = C 

is satisfied 

From the inelastic cross section 

v\n - 1 - (1 - oT{b)) A = £(oT(&))*(l - oT(b)) A - n (1) 

71=1 

retaining in the sum 

n 

(aT(b)r = Y,{ n t )°h<obT(bT (2) 

i=0 

the terms with i > 1, the C-cross section is obtained: 

a h c A = 1 - (1 - a c T(b)) A (3) 

If ac is small, <t^. a behaves as Aac, otherwise it behaves as A 2 ^ac- 
Let us next consider the particle distribution associated to a rare event C 
(Drcll-Yan pairs, annihilation cross section, particles of a given type produced 
in a rapidity interval, slow nucleons, high pt particle, etc.) 

If ac is the probability of event C to occur in an elementary collision and 
N(u) is the number of events with v elementary collisions in a nucleus- nucleus 
collision, the number of events where event C occurs is given by 

Nc(v)=Y.{\) ai c^-*cT- t N{ V ) (4) 
»=i 

If event C is rare (ac small), we can retain the leading term in ac and 

N c {v) ^ a c vN{v) (5) 
The total number of C-events and total number of events are: 
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^N c {v) = a c {v)N 



(6) 



Therefore 

P C {V) = -j^P{ V ) (7) 

The multiplicity distribution on produced particles P(n) is given by 

p ( n ) = J2 Yl ^{v)p{ni)p{n 2 )...p{n v ) (8) 

u—l ni+...+u^-n 

and the different moments of the multiplicity distribution can be expressed in 
terms of the different moments of the distribution on the number of elementary 
collisions and of the distribution on the particles produced in the elementary 
collision. In particular: 



(n) = {u)n 
D 2 _ {u 2 ) - {vf 
{n) 2 (u) 2 (v) n 2 



j^yj increases with the complexity of the systems involved, from 0.09 in 

e+e~ to 0.25-0.3 in pp collisions and to 0.8-1 in AB collisions. In general, the 
main contribution to the different normalized moments of the total multiplicity 
distribution is given by the normalized moments of the distribution on the 
number of elementary collisions, in such a way that in (0) v can be exchanged 
by n or Et, the transverse energy (experimentally n depends linearly on Et)' 



Pc(n) = ^yP(n), P C (E T ) = j^P(Et) (10) 

Relations (^) are universal, independent of ac- In figure 1 our prediction 
r Drell-Yan pairs according to ( |Io|) is compared with the NA38, SU data 
Similar good agreements are obtained for many other cases like high px 
particles produced in aa collisions, and Z° production in pp collisions 
(see figure 2) or pp annihilations, as can be seen in reference 1. 
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Figure 1: NA38 experimental associated multiplicity to Drell-Yan pairs in S-U collisions 
(cross points) compared to < ^ > P( n ) where P(n) is the experimental multiplicity distribu- 
tion of S-U collisions (squared points). 



2.1 J/ty Suppression 

In the case of J/^ production, final state destructive absorption will eliminate 
the linearity dependence of Nc{v) on v by making the effective number of 
collisions where event C appears smaller. This means that instead of equation 
(||) we have 

N c {v)=a c v £ N{v) , e < 1 (II) 

and 

Pc(Et) = J]^P(Et) (12) 

In figure 3 it is seen that the Et distribution associated to J/^ in SU 
collisions can be described by formula (|lj) with e = 0.7. 
Notice that from (0) and (0) 

Nj^(E T )/N DY (E T ) ~ 1/E T ~ S (13) 

This ratio decreases with Et and the curvature, the second derivative, is 
always positive. In all computations of J/\f r absorption, including destruction 
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Figure 2: Prediction for the associated multiplicity distribution for and Z° events (round 
black points) together with the experimental data (W^ cross points and Z° white stars) 
and the minimum bias distribution (squares). 



by comovers, the tendency for a large Et saturation occurs, which implies 
positive curvature. 

There is, however, another possibility for changing the v linearity of equa- 
tion (||): in a Quark Gluon Plasma (QGP), the J/\P and formation will be 
prevented!! Now ac is a function of v vanishing for large v. If the transition 
is very sharp: 



aj/v(v)=0 v>v* (14) 
or making a more reasonable approximation 

aj/*M = aj/^exp(-^ 2 /^* 2 ) (15) 

Now 

Nj^(E T )/N DY (E T ) ~ exp(-£;2/£;* 2 ) ( 16 ) 

which presents an essential difference to (|l3|). For Et > E T the ratio has a 
negative curvature. The experimental ratio for SU does not present any change 
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Figure 3: Experimental J/\E> over DY pairs as a function of Et in S-U collisions. 

in the curvature. At first sight, the NA50 experimental data in presents 
a change in the curvature for central events. However, it has been shown that 
it can be described without any changeEl using formula ([l^) with e = 0.6. This 
value is obtained previously from the low Et associated probability. 

Coming back to our basic relation ( |l6| ) we can mention another possibility 
for distinguish absorption or comovers effects from QGP formation. In the first 
case Pj^{Et) > P(Et) when Et — > oo, while if plasma is formed Pj^,{Et) < 
P(Et). The SU data are consistent with the first inequality. 

Independently, whether or not, QGP is formed, it is worth to ask for some 
reason to distinguish central PbPb collisions al SPS energies from the rest of 
collisions. The answer is percolation of stringsH. 

3 Percolation of Strings 

In many models of hadronic collisions, color strings are exchanged between 
projectile and target. The number of strings grows with the energy and with 
the number of nucleons of the participating nuclei. When the density of strings 
becomes high, the string color fields begin to overlap and eventually individual 
string may form a new string which has a higher color charge at its ends 
corresponding to the ST/ (3) sum of the color charges located at the ends of the 
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Figure 4: Percolation of strings in transverse space. 



original strings. As a result, the higher cploiLmeans a larger string tension and 
heavy flavor is produced more efficiently EjIIj. Also, as the effective number of 
strings decreases, the mean multiplicity is less than in the case of independent 
string fragmentation. 

In impact parameter space these strings are seen as circles inside the to- 
tal interaction area. As the number of strings increases, more strings overlap. 
Above a critical density of strings percolation occurs, so thai-paths of over- 
lapping circles are formed through the whole collision areaBfca, see figure 4. 
Along these paths the medium behaves like a color conductor. The percolation 
is a second order phase transition and give rise to a non-thermalized QGP on 
a nuclear scale. The percolation threshold r\ c is related to the critical density 
of circles n c by the expression 

rj c = nr 2 n c (17) 

where irr 2 is the transverse surface of an string. 

T] c has been computed using Monte-Carlo simulation, direct-connectedness 
expansions and other methods. All the results are in the range rj c = 1.12 — 1. 18. 
Taking r = 0.2 fm which reproduces the A enhancement observed in AB 
collisions □, it is obtained 

n c = 8.9 - 9.3 strings//m 2 (18) 

In table 1 the number of strings exchanged for central pp, SS, FeAir, SU 
and PbPb collisions is shown together with their densities. It is seen that at 
SPS energies only the density reached in central PbPb collisions is above the 
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Table 1: Number of strings (upper numbers) and their densities (fm 2 ) (lower numbers) in 
central pp, SS, SU and PbPb collisions at SPS, RHIC and LHC energies. 



VS(AGeV) 


Collisions 




PP 


SS 


FeAir 


SU 


PbPb 


19.4 


4.2 


123 


89 


268 


1145 




1.3 


3.5 


4.42 


7.6 


9.5 


200 


7.2 


215 


144 


382 


1703 




1.6 


6.1 


7.16 


10.9 


14.4 


5500 


13.1 


380 


255 


645 


3071 




2.0 


10.9 


12.67 


18.3 


25.6 



critical density. At RHIC, for central AgAg collisions, the density of strings is 
9, just at the critical value. At LHC energies, SS central collisions are already 
over the critical density. Let us mention that for central FeAir collisions, 
the critical density is reached at LHC energies. In the laboratory frame, this 
means that between 10 17 eV and 10 18 eV the percolation of strings occurs in 
central FeAir collisions. Several experiments of cosmic rays have pointed out 
a change on the behavior of the atmospheric cosmic ray cascade around the 
energies mentioned above. The maximum depth of the cascade becomes larger 
indicating that the energy of the primary is dissipated slower. This is what it 
was expected if percolation occursEj. The multiplicity of FeAir is suppressed 
and therefore the depth of the cascade increases. The usual belief is that a 
change of the chemical composition of the primary occurs, and as the energy 
increases in the range 10 16 eV to 10 19 eV there are less Fe and more protons 
in the primary cosmic ray. Instead of this we pointed out the possibility that 
the hadronic phase transition takes place without requiring any ad-hoc change 
of the chemical composition of the cosmic ray. 

Let us point out that if percolation occurs, the energy-momentum of the 
large cluster of strings is the sum of the energy-momentum of each individ- 
ual string, and therefore when this object fragments it will produce particles 
with very large longitudinal momentum, outside of the nucleon-nucleon kine- 
matical limits. This cummulative and spectacular effect could be checked by 
Brahms detector at RHIC. The percolation in this way behaves like a powerful 
accelerator. 

In 2-dimensional percolation it is known that tha-fraction <f> of the total 
area occupied by strings is determined by the formulaE-3 

cj> = 1 - exp(-r?) (19) 
In the vicinity of the phase transition it is satisfied the scaling law for the 
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number of clusters with n strings 



(v n ) =n T F(n a {ri~ri c )) |r?-r? c |<l n>l (20) 

where the critical indices are t = 187/91 and a = 36/91. The function F(z) 
is finite at z = and falls off exponentially for \z\ — » oo. The equation ( ^o|) is 
of limited value, since near rj =_n c the bulk of the contribution is still supplied 
by low values of n, for which ( pO|) is not valid. However from (|J) one can finds 
non-analytical parts of other quantities of interest at the transition point. In 
particular, one finds a singular part of the total number of clusters M = ^v n 
as 

A(.M) = c\rj - ?7 C | 8 / 3 (21) 

This singularity is quite weak: not only (M) itself but also its two first 
derivatives in r] stay continuous at 7/ = r\ c and only the third blows up as 
1 77 — r/ c | -1 / 3 . This singularity implies also a change on the behavior of the mul- 
tiplicity. The strength of this singularity would depend on the dependences of 
the multiplicity of a cluster on the number of strings. Reasonable assumptions, 
imply that the second derivative of the multiplicity blows up. 

4 Conclusions 

It is well known that the structure of multiple scattering and the shadowing can 
explain many unexpected features of experimental data. Here, as an example 
of this, it is shown that the self-shadowing events have an universal associated 
multiplicity distribution independently of the detailed type of events, as far as 
their cross section is small. 

Although most of the experimental data of present and may be future 
heavy ion collisions can be explained by conventional physics, it is sure that 
there will be surprises and excitement related to the forthcoming RHIC and 
LHC. As an example, we discuss the possibility of percolation of strings. If 
this occurs, exotic phenomena should appear like the production of many very 
fast particle with longitudinal momentum much larger that the permitted by 
the nucleon-nucleon kinematics or a dumping on the expected multiplicities. 
It is pointed out, that the physics of cosmic-ray and the new accelerators 
begin to overlap. In this sense, the observed change in the depth of maximum 
atmospheric shower could be related to the percolation of strings formed in 
iron- air collisions. 

In conclusion we thank N. Armesto, M. Braun and J. Dias de Deus who 
participated in different parts of the work reported here. We would like to 
thank the organizers of this workshop for such a nice meeting and the CICYT 
of Spain for financial support. 
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